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ABSTRACT 

We report the discovery of Lyman-alpha emission (Lya) in the bright galaxy EGSY-2008532660 
(hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by 
Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric 
redshift ( z p hot. = 8.57+0.43), apparent brightness (Hi 6 o = 25.26 ± 0.09) and red Spitzer/IRAC [3.6]- 
[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total 
integration of ~4.3 hours, our data reveal an emission line at ~11776 A which we argue is likely Lya 
at a redshift z spec = 8.6831 q;qo 4, in good agreement with the photometric estimate. The line was 
detected independently on two nights using different slit orientations and its detection significance 
is ~ 7.5cr. An overlapping skyline contributes significantly to the uncertainty on the total line flux 
although the significance of the detected line is robust to a variety of skyline-masking procedures. 

By direct addition and a Gaussian fit, we estimate a 95% confidence range of 1.0-2.5x 10~ 17 erg s _1 
cm -2 , corresponding to a rest-frame equivalent width of 17-42 A. EGSY8p7 is the most distant galaxy 
confirmed spectroscopically to date, and the third luminous source in the EGS field beyond z p hot > 7.5 
with detectable Lya emission viewed at a time when the intergalactic medium is believed to be fairly 
neutral. Although the reionization process was probably patchy, we discuss whether luminous sources 
with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter 
population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright 
sources promises important insight into the early formation of galaxies. 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution — galaxies: 
formation 


1. INTRODUCTION 

Our understanding of the early Universe has improved 
considerably in recent years through the photometric dis¬ 
covery of large numbers of high-redshift galaxies in both 
deep and gravitationally-lensed fields observed with the 
Hubble Space Telescope (e.g. Ellis et al. 2013; McLure 
et al. 2013; Bradley et al. 2014; Oesch et al. 2014; Zheng 
et al. 2014; Bouwens et al. 2015b). Although uncertain¬ 
ties remain, the demographics and limited spectroscopic 
follow-up of this early population has been used to ar¬ 
gue that star-forming galaxies played a significant role in 
completing the reionization of the IGM by redshift z ~ 6 
(e.g. Robertson et al. 2013, 2015; Kuhlen & Faucher- 
Giguere 2012; Bouwens et al. 2012, 2015a; Finkelstein 
et al. 2012, 2015; Ishigaki et al. 2015). A key observa¬ 
tion delineating the end of the reionization epoch is the 
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marked decline beyond z ~6 in the visibility of Lyman- 
alpha (Lya) emission seen in continuum-selected galaxies 
(Stark et al. 2010; Fontana et al. 2010; Ono et al. 2012; 
Pentericci et al. 2011; Schenker et al. 2012, 2014). As 
Lya represents a resonant transition it is readily scat¬ 
tered by the presence of neutral gas and thus acts as a 
valuable proxy for the state of the IGM. 

There are, however, several limitations in using Lya 
as a probe of reionization. Firstly, converting the declin¬ 
ing visibility of Lya into a neutral gas fraction involves 
complex radiative transfer calculations and several un¬ 
certain assumptions. The currently observed decline (e.g. 
Schenker et al. 2014) implies a surprisingly rapid end 
to the process (Choudhury et al. 2014; Mesinger et al. 
2015). Secondly, simulations suggest that reionization 
is likely to be a patchy process (Taylor & Lidz 2014) 
and thus conclusions drawn from the modest samples 
of spectroscopically-targeted galaxies may be misleading 
(Treu et al. 2012; Pentericci et al. 2014; Tilvi et al. 2014). 
Finally, the absence of a Lya detection in the spectrum 
of a proposed high-redshift candidate may simply imply 
the source is a foreground galaxy (although the contam¬ 
ination rate for such galaxies is typically very low, e.g. 
Vanzella et al. 2009; Stark et al. 2010). Ideally targets 
in such studies would be spectroscopically confirmed in¬ 
dependently of Lya, for example using UV metal lines 
(Stark et al. 2015b, a, see also Zitrin et al. 2015). 

Although observational progress is challenged by the 
faintness of targets selected in deep fields such as the 
Hubble Ultra Deep Field (typically wiab ~27-28), an im- 
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Fig. 1.— Spectroscopic detection of emission in EGSY8p7 with MOSFIRE. Upper panel shows the 2D spectrum below which we plot 
the raw ( black line) and smoothed ( blue line) ID spectrum and its error ( red shading). The red line shows an example best-fit model of 
the data (§3). Vertical lines mark OH skyline positions. The upperleft panel shows a normalized signal map extracted along the slit within 
a 5-pixel (~ 6.5 A) wide box centered on the line. The pattern of two negative peaks bracketing the positive peak exactly matches that 
expected from the dithering scheme used. Arrows show the predicted locations of other lines for a lower redshift interpretation of the line. 
Green boxes on the 2D spectrum mark the skyline region typically masked out in our calculations. See §3 for more details. 


portant development has been the identification of much 
brighter z > 7 candidates from the wider area, some¬ 
what shallower, Cosmic Assembly Near-infrared Deep 
Extragalactic Legacy Survey (CANDELS, Grogin et al. 
2011; Koekemoer et al. 2011). Surprisingly, some of these 
brighter targets reveal Lya despite lying inside the pu¬ 
tative partially neutral era. Finkelstein et al. (2013, 
hereafter F13) reported Lya with a rest-frame equiv¬ 
alent width (EW) of 8 A at £=7.508 in a 77,45=25.6 
galaxy; Oesch et al. (2015, hereafter 015) find Lya emis¬ 
sion at £=7.73 with EW=21 A in an even brighter source 
at 77,45=25.03; and Roberts-Borsani et al. (2015, here¬ 
after RB15) identified a tentative Lya emission (4.7er) 
in a 77,45=25.12 galaxy at a redshift £=7.477, which we 
have now confirmed (Stark et al, in prep). In addition 
to their extreme luminosities (Mjjv — — 22), these three 
sources have red [3.6]-[4.5] Spitzer/IRAC colors, indica¬ 
tive of contamination from strong [O III] and Balmer H/3 
emission. 

Using the Multi-Object Spectrometer For Infra-Red 
Exploration (MOSFIRE, McLean et al. 2012) on the 
Keck 1 telescope, we report the detection of a promi¬ 
nent emission line in a further bright candidate drawn 
from the CANDELS program. EGSY-2008532660 (here¬ 
after EGSY8p7; RA=14:20:08.50, DEC=+52:53:26.60) 
is a 77,45=25.26 galaxy with a photometric redshift of 
8 .57lg;43 an d a re d IRAC [3.6]-[4.5] color, recently dis¬ 
covered by RB15. We discuss the likelihood that the line 


is Lya at a redshift z spec = 8.68 making this the most dis¬ 
tant spectroscopically-confirmed galaxy. Detectable Lya 
emission at a redshift well beyond £ ~ 8 raises several 
questions regarding both the validity of earlier claims for 
non-detections of Lya in fainter sources, and the physical 
nature of the luminous sources now being verified spec¬ 
troscopically. Even if these bright systems are not rep¬ 
resentative of the fainter population that dominate the 
ionization budget, they offer new opportunities to make 
spectroscopic progress in understanding early galaxy for¬ 
mation. 

The paper is organized as follows: In §2 we review 
the object selection, spectroscopic observations, and data 
reduction. The significance of the line detection and its 
interpretation as Lya is discussed in §3. We discuss the 
implications of the detectability of Lya in the context of 
the earlier work in §4. Throughout we use a standard 
ACDM cosmology with f l m0 = 0.3, SIao = 0.7, 77 0 = 100 
h km s -1 Mpc , h = 0.7, and magnitudes are given 
using the AB convention. Errors are lcr unless otherwise 
stated. 

2. DATA 

The galaxy EGSY8p7 was detected in the Extended 
Groth Strip (EGS; Davis et al. 2007) from deep (> 27.0) 
multi-band images in the CANDELS survey and first re¬ 
ported as one of four unusually bright (77 i 60 < 25.5) 
candidate z > 7 galaxies by RB15. One of these, EGS- 
zs8-l, with Zphot = 7.92±0.36 was spectroscopically con- 
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Night 1 




Night 2 



Fig. 2.— Confirmation of line detection in EGSY8p7 over two nights. Left: HST F160W image of the EGSY8p7 field with the slit 
orientations adopted in two successive nights. Right: Extracted ID signal/noise spectra across the J-band for each night, with a zoom of 
the 2D data around the line, marked with a vertical blue line. The 2D spectra are smoothed with a 3-pixel Gaussian for better illustration 
of the data. The ID spectrum is smoothed with a Gaussian of a = 5 A, comparable to the measured line width. The Y axis is scaled 
so that the peak signal/noise matches the integrated value (e.g. F13). Horizontal dashed lines mark the ±3(7 region. On both nights, the 
signal/noise at the line location significantly exceeds that elsewhere. 


firmed at z = 7.73 by 015. Normally such objects would 
be selected as Y-band dropouts but, given Y-band ob¬ 
servations are not yet available over the full CANDELS 
field, an alternative selection criterion was adopted that 
takes advantage of red IRAC [3.6]-[4.5] colors, indicating 
prominent [O IIIJ+H/3 emission within the 4.5 /im band 
(see also Labbe et al. 2013; Smit et al. 2015). Of the ob¬ 
jects listed by RB15, EGSY8p7 is the only one for which 
Y-band data is not available. As a result, its photometric 
redshift derived from HST data alone is fairly uncertain 
(5.6 < z p hot < 9.2; Fig. 5 of RB15) but including its 
IRAC [3.6]-[4.5] color of 0.76±0.14 (and blue JH 140 -[3.6] 
color) narrows the range to z p hot = 8.57^43 ■ 

We observed EGSY8p7, the highest-redshift candidate 
in RB15’s list, on June 10 and 11 2015 with MOSFIRE 
on the Keck 1 telescope. Observations in the J-band 


spanned the wavelength range 11530A < A < 13520A 
using an AB dithering pattern of ±1.25" along the slit 
with individual frames of 120s. The slit masks on the 
two nights differed by 120° in orientation both with a 
slit width of 0.7". In each mask we allocated one slit 
to a nearby star to monitor changes in seeing, trans¬ 
parency, and possible positional drifts. Conditions were 
clear throughout, with an average seeing of 0.60" for the 
first night and 0.76" for the second night. The total 
nightly exposure times were 158 minutes and 128 min¬ 
utes. Excluding exposures where the seeing was signif¬ 
icantly (> la ) worse than average, the useful exposure 
time comprises all frames from the first night and 80% of 
those acquired on the second night (where the conditions 
were less stable) culminating in a total exposure-time of 
4.33 hours. Calibrations were obtained via long-slit ob- 
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servations of standard AOV stars. 

Data reduction was performed using the standard 
MOSFIRE reduction pipeline 9 . For each flat-Eelded slit 
we extracted the ID spectrum using a 11-pixel boxcar 
centered on the expected position of the target. A sim¬ 
ilar procedure was adopted in quadrature to derive the 
lfj error distribution. To ensure that the derived er¬ 
ror spectrum reflects the noise properties of the data, 
we measured the standard deviation of the pixel-by-pixel 
SNR, which should be unity for a set of independent val¬ 
ues drawn from Gaussian distributions. We obtained a 
value of 1.3, implying that the errors are slightly under¬ 
estimated by the pipeline, and corrected the error spec¬ 
trum accordingly. Data from both nights were co-added 
by inverse-variance averaging the calibrated ID spectra. 
We also allowed for relative shifts along the slit of 2 pix¬ 
els (~0.2") and across the slit of 0.2"(which affects the 
expected slit loss correction) and propagated the associ¬ 
ated uncertainty into our error budget. All reductions 
and calibration steps were performed independently by 
two authors (AZ, SB). 

To calibrate the spectra we scaled a Vega model 10 to 
each standard star to determine a wavelength-dependent 
relative flux calibration and telluric correction using the 
procedure described in Vacca et al. (2003). Independent 
telluric calibrations derived from three standard stars 
agree to within 5% at the location of the detected emis¬ 
sion line. The absolute calibration for each night was de¬ 
rived by comparing the spectroscopic magnitude of the 
stars on the slitmask with photometric measures in the 
J-band obtained from the 3D-HST catalog (Brammer 
et al. 2012; Skelton et al. 2014). Note that this procedure 
takes into account differential slit losses due to varying 
seeing. The absolute calibration factors obtained for the 
two nights are consistent with differences arising from 
seeing effects. We adopt a conservative absolute calibra¬ 
tion error of 30%, taking into account the nightly vari¬ 
ation, differences between the independent reductions, 
and other contributions mentioned above. The final 2D 
and ID spectra are shown in Fig. 1. 

3. LYMAN a AT A REDSHIFT 8.68 

Fig. 1 reveals a prominent line at A ~ 11776 A flanked 
by a skyline on its blue side. Examining the signal and 
associated noise in a circular aperture of 6 pixels in radius 

- corresponding to ~ 2 times the line width found below 

- we estimate a significance of 7.6cr in the 2D spectrum 
and, within the extracted ID spectrum over the same 
spectral range, 7.5cr. This significance holds over a range 
of integration wavelengths both including and excluding 
the skyline. Additionally, the line is detected indepen¬ 
dently on each night (6.0cr on night 1, 3.8<r on night 2), 
despite the changed slit orientation (Fig. 2). Taking the 
signal/noise on the first night, and assuming it scales 
as Vf/FWHM, where t is the exposure and FWHM the 
seeing, for the second night (where the exposure was less 
and seeing worse) we predict a signal/noise similar to 
that observed. 

We estimate the observed line properties by using a 
Monte Carlo Markov Chain to fit a truncated Gaussian 

9 http://www2.keck.hawaii.edu/inst/mosfire/drp.html 

10 http://kurucz.harvard.edu/stars/vega/ 


TABLE 1 

Emission Line properties 


Name Fiducial model 95% C.I. 


f(Lya) a 

1.7 

[1.0 - 2.5] 

/A 

11774 A 

[11770 - 11776] A 

z c 

8.683 

[8.679 - 8.684] 

EW d 

28 A 

[17 - 42] A 

a e 

4.7 A 

[2.0 - 8.4] A 

Vj 

118 km s -1 

[51 -214] km s” 1 

Vfwhm s 

277 km s” 1 

[120 -503] km s" 1 


Note. — 

Emission line properties derived from a truncated Gaus¬ 
sian fit (Fig. 1), corrected for instrumental broadening, 
and not corrected for IGM absorption. The 95% range 
quoted is adopted from a set of fits as described in the 
text. 

a Total flux in units 10 -17 erg s^ 1 cm^ 2 
° Peak of Gaussian 
c Redshift 
d Rest-frame EW 
e Gaussian line width 

f Velocity width, corresponding to the Gaussian width 
g Velocity FWHM (~ 2.35 X a) 

to the data (in a similar way to 015), taking into ac¬ 
count instrumental broadening. We leave the truncation 
point a free parameter, and mask out pixels contami¬ 
nated by strong skyline residuals. The proximate skyline 
significantly affects the fit, and different masking con¬ 
figurations yield different results. To span the range of 
possible solutions we repeat the model fitting, as well as 
direct integration measurements, with a variety of mask¬ 
ing configurations. In each such configuration a differ¬ 
ent part of the skyline was excluded. For each of these 
resulting models we calculate the 95% confidence level 
intervals. For each parameter we then take the union of 
the intervals from different models and report them in 
Table 1. 

The line may be asymmetric with some attenuation on 
its blue side, but this is difficult to ascertain given the 
proximate skyline, and the best-fit models (Fig. 1) are 
mostly symmetric. The line peak lies between 11770 and 
11776 A which, for Lya, implies a redshift of 8.683^°'oo4- 
Note the typical peak wavelength is a few A bluer than 
the raw peak in Fig. 1. The line width ranges from 2.0- 
8.4 A, corresponding to a velocity broadening, corrected 
for instrumental effects, of V CT =51-214 km s . The line 
flux ranges from 1.0-2.5xl0 -1 ' erg s _1 cm -2 (excluding 
the 30% uncertainties in the absolute calibration). Im¬ 
portantly, the significance of the line remains at the > 7er 
level for all chosen integration methods and variations in 
the masking of the skyline. The reliability of the SNR 
measurement given the nearby skyline was further veri¬ 
fied using simulated lines inserted into the 2D spectra. 

We can estimate the chance of finding such a line by 
examining the signal/noise distribution in our combined 
2D spectral data at 1500 random locations taking the 
same aperture used for measuring the line properties 
above. The resulting histogram in Fig. 3 reveals a very 
low probability of finding such a feature by chance.Our 
measured signal/noise represents a 6er deviation in the 
distribution (which has a width larger than unity due 
to correlated noise). We also verified that proximity to 
skylines does not increase the width of the distribution. 
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Fig. 3. — Verifying the reality of the line detection in EGS8p7. 
The signal/noise distribution from 1500 randomly-positioned aper¬ 
tures across the 2D data presented in Fig. 1. All apertures were 
6 pixels in radius (~ 7.5 A in wavelength space), as the aperture 
used to measure the line significance. The vertical line shows the 
detected signal/noise far exceeds that of the noise distribution. 

Together with the facts that the positive and negative 
line positions in the 2D spectrum match the dithering 
pattern used (Fig.l) and the line position is astromet- 
rically centered on the target position along the slit on 
both nights, we conclude the probability for the detection 
to be an artifact is negligibly small. 

Adopting the identification of Lya, the inferred red- 
shift of z = 8.683lg;oo4 is satisfactorily close to the pho¬ 
tometric estimate, z p hot = 8.57+0(43 reported by RB15. 
Using the HST iA-band photometry unaffected by line 
emission, for a reasonable UV slope (/3 = —2), the range 
of line fluxes (Table 1) translates to a rest-frame EW of 
17-42 A. While formally fairly uncertain, it is compa¬ 
rable to line properties secured with MOSFIRE for the 
other bright galaxies at z p hot > 7.5 selected in CAN- 
DELS with red IRAC colors (F13, 015). 015 obtained 
a total line flux of 1.7 ± 0.3 x 10 -17 erg s -1 cm -2 and 
a rest-frame EW of 21 ± 4 A, whereas F13 secured a to¬ 
tal line flux of 2.64 x 10 -18 erg s _1 cm -2 and a more 
modest rest-frame EW~7.5 A. We note RB15 discussed 
the possibility that EGSY8p7 is magnified by foreground 
galaxies. Adopting the lens galaxy masses found in RBI5 
as input, and using an updated photometric redshift, we 
conclude EGSY8p7 is likely magnified by ~ 20%, and 
less than a factor ~ 2. 

Given we only detected one spectral line in the J-band, 
we also considered interpretations other than Lya, corre¬ 
sponding to a lower redshift galaxy. If the line is a com¬ 
ponent of the [O II] (AA3726, 3729) doublet, this would 
imply z ~ 2.16, possibly consistent with a low redshift 
solution to the SED (Fig. 4; z p hot = 1.7±0.3; 99% C.I.). 
However, for each [O II] component we find no trace of 
the other line to a limit several times fainter than the de¬ 
tected line ( magenta arrows in Fig. 1). Similarly, in the 
case of [O III] (AA4959, 5007; cyan arrows), or the Hg 
line (A4861), implying z ~ 1.4, we would expect to see 
(at least) one of the other lines given typical line ratios 
assumed. While the presence of skylines clearly limits a 



Fig. 4.— Photometric SED fitting strongly supports a high- 
redshift solution for EGSY8p7: Red points show magnitudes, 
where the horizontal error marks the filter width, and triangles 
denote la upper-limits from non-detections. The blue SED shows 
the best-fit template at 2=8.57. The magenta SED shows the 
z=8.48 solution adopting the strong emission line template used in 
RB15 (see their work for more details on the templates). The grey 
SED shows the best-fit low redshift template which cannot explain 
the absence of detection in the optical bands and which is ~ 10 7 
times less likely than the high-redshift solutions. 

robust rejection of these alternatives, both the absence 
of any optical detections and the red IRAC color are 
hard to understand in a low redshift galaxy with intense 
line emission. Formally, the SED fit, performed using 
the EAZY program (Brammer et al. 2008), gives a low 
redshift likelihood 10 7 times smaller than our adopted 
solution at 2 = 8.68 (Fig. 4). 

4. DISCUSSION AND SUMMARY 

Although we claim EGSY8p7 is the highest-redshift 
spectroscopically-confirmed galaxy and the first star¬ 
forming galaxy secured beyond a redshift of 8 (c.f. a 
gamma ray burst at z = 8.2, Tanvir et al. 2009), 
its prominent Lya emission, detected in only 4 hours, 
raises several interesting questions about the usefulness 
of galaxies as tracers of reionization. The last few years 
has seen a consistent picture emerging from observational 
programs charting the visibility of Lya beyond a red¬ 
shift z ~6.5 in much fainter (m ~ 27) galaxies (Schenker 
et al. 2014 and references therein). Although the expo¬ 
sure times, selection techniques and observational strate¬ 
gies differ across the various programs, almost 100 faint 
targets have now been targeted, collectively yielding very 
few z > 7 Lya detections (F13, Schenker et al. 2014). In 
marked contrast, all 3 bright z p hot <1 7.5 targets in the 
EGS field in Table 2 of RB15’s compilation now reveal 
Lya emission. 

Although this bright spectroscopic sample is admit¬ 
tedly modest, the high success rate in finding Lya at 
redshifts where the IGM is expected to be substantially 
neutral (~70% according to Robertson et al. 2015) raises 
two interesting questions. Firstly, do these new results 
challenge claims made by the earlier fainter Lya searches 
in the same redshift interval? Is it possible that those 
campaigns reached to insufficient depths to reliably de- 
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tect Lya emission? Typically such studies adopted the 
rest-frame EW distribution at z ~6 published by Stark 
et al. (2011) and simulated, using their appropriate sur¬ 
vey and instrumental parameters, the likely success rate 
to somewhat higher EW limits of ~25-50 A. A key issue 
in understanding whether the RB15 objects are similarly 
attenuated by the partially neutral IGM at z > 7.5 is 
the fiducial z ~6 EW distribution for these luminosities. 
Stark et al’s EW distribution is based on an analysis of 
74 2 ~ 6 galaxies targeted spectroscopically within the 
deep GOODS fields but only 13 reached the UV luminosi¬ 
ties being probed here. Curtis-Lake et al. (2012) present 
a smaller sample of 10 spectroscopically-confirmed z ~6 
galaxies from the UDS field, several of which are more 
luminous with EWs > 40 A. A key deficiency in this 
comparison is that the z ~ 6 samples were not selected 
to have as large EW [OIIIJ+H/3 emission implied in the 
IRAC [3.6]-[4.5] color-selected RB15 sample. Stark et al. 
(2014) show that Lya EWs are uniformly extremely large 
(40 -160 A) in 2 ~2-3 galaxies with such large EW op¬ 
tical line emission. Thus it is possible that Lya in these 
Zphot > 7.5 galaxies is attenuated. In addition to the 
possibility of cosmic variance in the EGS field, it may be 
that such luminous objects at 2 > 7 lie in large overdense 
regions and present a somewhat accelerated view of the 
reionization process (Barkana & Loeb 2006). 

The second question is whether these rarer luminous 
sources (c.f. Matthee et al. 2015; Ouchi et al. 2013) are 
physically distinct from the fainter population. Adopt¬ 
ing the IRAC red [3.6]-[4.5] color selection may preferen¬ 
tially select unusually intense line emission with abnor¬ 
mal ionizing spectra capable of creating early bubbles of 
ionized gas in the local IGM (see also F13, Ono et al. 
2012). Evidence for stronger than usual ionizing spectra 
in 2 > 7 dropouts with known Lya emission is provided 
by the prominent detection of CIV 1548 A emission in 
the lensed galaxy at 2 = 7.047 (Stark et al. 2015a). The 
key issue in this case is the proportion of 2 > 6.5 Ly¬ 
man break galaxies that reveal such prominent IRAC 
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excesses. Smit et al. (2015) demonstrate that perhaps 
upwards of 50% of the photometrically selected galaxies 
at z ~6.6-6.9 have IRAC excesses that are comparable 
in magnitude to those seen in RB15. 

In summary, not only have we pushed the detectability 
of Lya emission to well beyond a redshift z ~8, but it 
seems luminous z p i lot > 7.5 galaxies selected additionally 
for their red IRAC [3.6]-[4.5] colors have an unusually 
high success rate of line detection. We present several 
reasons why this might contrast with the low yield of de¬ 
tecting Lya in intrinsically fainter sources. Regardless 
of the explanation, it is clear that further spectroscopic 
follow-up of such examples will yield valuable informa¬ 
tion on the early development of massive galaxies. 

ACKNOWLEDGEMENTS 

We thank the reviewer of this work for useful com¬ 
ments. We acknowledge useful discussions with Jim Dun¬ 
lop, Brant Robertson and Chuck Steidel. We thank 
Garth Illingworth, Pieter van Dokkum, Iva Momcheva, 
Marijn Franx, Mauro Stefanon, and Benne Holwerda for 
their role in the selection of this object. Support for this 
work was provided by NASA through Hubble Fellowship 
grant #HST-HF2-51334.001-A awarded by STScI, which 
is operated by the Association of Universities for Re¬ 
search in Astronomy, Inc. under NASA contract NAS 5- 
26555. RSE acknowledges the hospitality of the Insti¬ 
tute for Astronomy, Edinburgh via a Carnegie Centen¬ 
nial Professorship. This work is in part based on pre¬ 
vious observations made with the NASA/ESA Hubble 
Space Telescope, the Spitzer space telescope, the Sub¬ 
aru telescope, and the Canada France Hawaii Telescope. 
The data presented herein were obtained at the W.M. 
Keck Observatory. The authors wish to recognize and ac¬ 
knowledge the very significant cultural role and reverence 
that the summit of Mauna Kea has always had within the 
indigenous Hawaiian community. We are most fortunate 
to have the opportunity to conduct observations from 
this mountain. 


Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 
197, 35 

Ishigaki, M., Kawamata, R., Ouchi, M., et al. 2015, ApJ, 799, 12 
Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, 
ApJS, 197, 36 

Kuhlen, M., Sz Faucher-Giguere, C.-A. 2012, MNRAS, 423, 862 
Labbe, I., Oesch, P. A., Bouwens, R. J., et al. 2013, ApJ, 777, L19 
Matthee, J., Sobral, D., Santos, S., et al. 2015, MNRAS, 451, 4919 
McLean, I. S., Steidel, C. C., Epps, H. W., et al. 2012, in Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 8446, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, 0 
McLure, R. J., Dunlop, J. S., Bowler, R. A. A., et al. 2013, 
MNRAS, 432, 2696 

Mesinger, A., Aykutalp, A., Vanzella, E., et al. 2015, MNRAS, 
446, 566 

Oesch, P. A., Bouwens, R. J., Illingworth, G. D., et al. 2014, ApJ, 
786, 108 

Oesch, P. A., van Dokkum, P. G., Illingworth, G. D., et al. 2015, 
ApJ, 804, L30 

Ono, Y., Ouchi, M., Mobasher, B., et al. 2012, ApJ, 744, 83 
Ouchi, M., Ellis, R., Ono, Y., et al. 2013, ApJ, 778, 102 
Pentericci, L., Fontana, A., Vanzella, E., et al. 2011, ApJ, 743, 132 
Pentericci, L., Vanzella, E., Fontana, A., et al. 2014, ApJ, 793, 113 
Roberts-Borsani, G. W., Bouwens, R. J., Oesch, P. A., et al. 

2015, arXiv, 1506.00854 


Spectroscopic Measurement of a Redshift z=8.68 Galaxy 


7 


Robertson, B. E., Ellis, R. S., Furlanetto, S. R., &; Dunlop, J. S. 
2015, ApJ, 802, L19 

Robertson, B. E., Furlanetto, S. R., Schneider, E., et al. 2013, 
ApJ, 768, 71 

Schenker, M. A., Ellis, R. S., Konidaris, N. P., & Stark, D. P. 
2014, ApJ, 795, 20 

Schenker, M. A., Stark, D. P., Ellis, R. S., et al. 2012, ApJ, 744, 
179 

Skelton, R. E., Whitaker, K. E., Momcheva, I. G., et al. 2014, 
ApJS, 214, 24 

Smit, R., Bouwens, R. J., Franx, M., et al. 2015, ApJ, 801, 122 
Stark, D. P., Ellis, R. S., Chiu, K., Ouchi, M., & Bunker, A. 2010, 
MNRAS, 408, 1628 

Stark, D. P., Ellis, R. S., &; Ouchi, M. 2011, ApJ, 728, L2 
Stark, D. P., Richard, J., Siana, B., et al. 2014, MNRAS, 445, 
3200 

Stark, D. P., Walth, G., Chariot, S., et al. 2015a, arXiv, 
1504.06881 


Stark, D. P., Richard, J., Chariot, S., et al. 2015b, MNRAS, 450, 
1846 

Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, Nature, 461, 
1254 

Taylor, J., & Lidz, A. 2014, MNRAS, 437, 2542 

Tilvi, V., Papovich, C., Finkelstein, S. L., et al. 2014, ApJ, 794, 5 

Treu, T., Trenti, M., Stiavelli, M., Auger, M. W., & Bradley, 

L. D. 2012, ApJ, 747, 27 

Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 
389 

Vanzella, E., Giavalisco, M., Dickinson, M., et al. 2009, ApJ, 695, 
1163 

Zheng, W., Shu, X., Moustakas, J., et al. 2014, ApJ, 795, 93 
Zitrin, A., Ellis, R. S., Belli, S., &; Stark, D. P. 2015, ApJ, 805, L7 



